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Abstract. The high-frequency-peaked BL Lac, MS 0205.7-1-3509 was observed twice with XMM-Newton. Both X-ray spectra 
are synchrotron-dominated, with mean 0.2-10 keV fluxes of 2.80 + O.OI and 3.34 + 0.02 x 10"'^ erg cm"^ s"' . The X-ray spectra 
are well fit by a power-law with absorption above the Galactic value, however no absorption edges are detected, implying a 
low metallicity absorber (Zq = 0.04;^q[5]) or an absorber with redshift above one (best-fit z = 2.1 for an absorber with solar 
abundances). In either case the absorbing column density must be ~ 9 x 10^' cm"^. A new optical spectrum is presented, with 
a Mg II absorption doublet detected at z = 0.351, but no other significant features. The optical spectrum shows little reddening, 
implying a low dust to gas ratio in the absorber MS 0205.7+3509 must therefore be viewed through a high column density, 
low-metallicity gas cloud, probably at z = 0.351 and associated with the galaxy that has been shown to be within ~ 2" of the 
BL Lac. 
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d 1. Introduction that BL Lacs are gravitationall y microlensed qu asars can be 

discounted for most BL Lacs ( |Padovanilll99a and only a 

Blazars are divided into BL Lacs and quasars (either flat- ^ • i-i u ui j j 4. . . ui 

^ ^ tew remain as possible or probable candidates, most notably 

spectrum radio-loud, optically violently variable, highly- »--,rmc- u * u f j i. 

^ ' t- J J ^ b J AO0235-H164 which appears to have foreground absorption 

polarised or core-dominated) based on the strength of emission \^^,^^^^ ^,^^\ \^ |Sd^ 1994, Madeiski^ S 

lines in the optical spectrum (Scarpa & Falomo 1997, and ref- , . a<^m ,n u j . 1 inn^^ r\^u j- 

. ^ ^ and a companion AGN (Burbidg e et al.. .l996). Other candi- 

erences therein). dates include PKS 0537-441 (which shows rapid microvariabil- 

It has been suggested that sources with BL Lac character- -^^^ ^^^^ ^^^^ evidenc e for a foregro und obje ct in op- 

istics are actuaUy gravitationaUy microlensed quasars ( Nottale tj^^j ^^^^^^ spectroscopy, iRomero et ab.l999. .PkT^ 

| l986l | Ostriker&Vietri|l99(]|) . In these cases, it is expected that ^q2), and B2 1308+326 which has c haracteristics in terme- 

stellar mass lenses in a foreground galaxy significantly ampli- ^iate betwee n BL Lacs and quasars JOabuzda et alJ fl99i 

fies the central QSO continuum source but not the emission |watson et aljEoOol) . but where high resolution imaging of the 

from the line-emitting regions and that variations in the rela- Lac with flie HST WFPC2 ( Urrvet al. 1999) was consistent 

five source-lens position could account for the rapid variabil- ^j^j^ ^ point source 

ity observed in many BL Lacs. Sources of this kind should 

clearly have foreground galaxies, which would result in an ap- 0205.7+3509 is another such rare candidate, and while 

parent decentering of the AGN from the "host" and an ex- ^^^^P ^^§^"8 ^as revealed that the BL Lac is centred on a host 

cess of absorption in these sources. However, the suggestion g^l^^^^ ^hat is likely an elliptical and not offset in a spiral host 

as had been inferred from previous observations (Sto cke et alJ 

Send ojfprint requests to: D. Watson, darach@astro .ku.dk F995^ . a companion galaxy was also detected very close to the 
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BL Lac line of sight jFalomo et al J 19971) which had caused the 
previous inference of decenterin g in a spiral host t o be made. 
X-ray observa tions with ROSAT JStocke et al.ll 995) wdASCA 
dWatson et alji,199 9) showed the existence of absorption well 
above the Galactic level, inde ed at a level second only among 
BL Lacs to PKS 1413+135 jPerlman etaf]l2002h . These re- 
sults, on MS 0205.7+3509, imply that the X-ray absorber is 
in the companion galaxy which is foreground to the AGN. It 
has been suggested that stars in the halo of this companion 
galaxy c ould be responsible for microl ensing of the BL Lac 
(Falomo etalJl997llwatson et al.ll999l) . 

In spite of the relatively good spectral resolution of ASCA, 
the redshift of the absorber could not be constrained from those 
observations ( Watson et al. 1999). A redshift of z = 0.318 
was proposed based on the tentative detection of a Can ab- 
sorption system reported in the optical spectrum of this source 
(Morris et al. 1991) indicating the possible redshift either of 
the host galaxy or of a foreground absorb er and to date this 
redshift had been used as the best available jRector et alJ2000l 
[Stocke et al. 1995). 

Though somewhat lower in terms of column den- 
sity, MS 0205.7+3509 is six times brighter in X-rays than 
PKS 1413+135, making it one of the best available cases in 
which to study absorption in the hot phase of the ISM in a 
galaxy that is not at low redshift. XMM-Newton observations 
were performed in an attempt to determine the nature of the X- 
ray absorber in MS 0205.7+3509, in particular in the context of 
a foreground lensing galaxy. 

Results from these data are presented in this paper Sect. 2 
deals with the observations and the data reduction procedures; 
results from the X-ray and optical data are in Sect. 3. A dis- 
cussion of these results and a summary of our conclusions 
are given in Sect. 4. Uncertainties given are 90% confidence 
limits unless otherwise stated. A flat universe with Hq - 
75kms ' Mpc ' and Qa = 0.7 are assumed throughout. 

2. Observations and data reduction 

MS 205.7+3509 was observed by XMM-Newton jjansen et all 
I2OOII) for40ks and 20 ks during orbits 217 and 395 with obser- 
vations starting on 14 Feb. 2001, 07:06:41 UT and 4 Feb. 2002, 
19:41 :06UT respectively. 

For the first observation, exposures of 38 ks and 34 ks d ura- 
tion were made with each EPIC-MOS ("Turner etal."2001^ and 
with the EPIC-pn (Striider et al. 2001) detector respectively; 
while for the second observation, the exposures were 17ks 
(each EPIC-MOS) and 15 ks (EPIC-pn). All observations were 
performed in Full Frame mode, using the Medium filter for the 
MOS cameras and the Thin filter for the pn. Data were pro- 
cessed and screened in a standard way with the XMM-Newton 
Science Analysis System version 5.3.0 (SAS). Only events cor- 
responding to patterns 0-12 were used for the two MOS cam- 
eras, while pattern 0-4 events were selected from the pn data. 

A circular extraction region was defined around the cen- 
troid position of the source, with an aperture of 35" radius. 
Data from these spatial regions were used to extract spectra 
and lightcurves. An estimate of the background was derived by 
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Fig. 2. Optical spectrum from the MMTO. Absorption lines of 
Fe II and Mg ii at z = 0.351 are labelled. Inset. A blow-up cen- 
tred on the resolved Mg ii doublet. Wavelength (A) is plotted on 
the axis of abscissas, flux (arbitrary units) on the ordinate. (The 
deviations from a power-law shape at the blue and red ends of 
the spectrum are due to calibration and atmospheric effects and 
affect only the continuum shape.) 

using an aperture of 100" radius at a source-free position close 
to the source extraction region. 

The complete lightcurves were created by binning the data 
in 70s bins (giving ~ 150 counts per bin); two other, different 
time binnings yielded similar results. The background for each 
instrument was subtracted from the source lightcurve and the 
per-instrument lightcurves were added. 

The extracted spectra were binned to give a minimum of 
twenty counts per bin. Data in the energy range 0.2-lO.OkeV 
were used. The background-subtracted spectrum from each in- 
strument was fit separately using the Levenberg-Marquardt al- 
gorithm in Xspec vl 1.0.1 and EPIC response matrices gener- 
ated with the SAS. 

The EPIC-pn 0.2-lO.OkeV count rate was found to be be- 
tween 1.2-1.5 counts/s; effects due to photon pile-up are negli- 
gible at this flux level (< 0.2% of the total counts). 

Differences between the spectra extracted from the EPIC- 
MOS cameras are negligible in this observation. The MOS data 
were therefore co-added and fit as a single spectrum, using the 
co-added response matrices and backgrounds. However, known 
differences in the cross-calibration of the EPIC-pn and MOS 
instruments are apparent in these spectra (Fig. these dif- 
ferences primarily affect the low-energy power-law slopes and 
normalisations. In order to exploit the full dataset, the MOS and 
pn spectra have been fit separately and then simultaneously and 
any significant differences have been highlighted. 

2.1. Optical 

An optical spectrum (Fig.|2j was obtained in a 9200 s exposure 
completed on 23 November 2000 with the Blue Channel spec- 
trograph at the MMTO.' The slit width was 1".25, resulting 

' The MMT Observatory is a joint facility of the University of 
Arizona and the Smithsonian Institution 
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channel enengy (keV) 



Fig. 1. Simultaneous fit 
: of EPIC-pn (crosses) and 

- combined EPIC-MOS (tri- 
angles) data from 14 Feb. 

! 2001 folded through the 

- detector response with fit 
residuals. The model fit 

E was a Galactic-absorbed 
; power-law with neutral 
1^ low-abundance absorption 
at redshift z = 0.351. The 
elemental abundance and 
_ absorbing column were 
fit simultaneously in both 
datasets (from 14 Feb. 
2001 and 4 Feb. 2002), 
whereas the power-law 

- slopes and normaUsations 
were allowed to vary in- 
dependently to account for 
the temporal variability of 
the synchrotron spectrum. 
Inset. Power-law model fit 
to the 1-lOkeV EPIC-pn 
data of 14 Feb. 2001 with 
Galactic absorption only. 



in a FWHM of 4.2 A. The slit position was E-W. The signal to 
noise ratio is ~ 20. The spectrum was extracted using the IRAF 
kpnoslit package. The spectra were wavelength calibrated with 
HeNeAr spectra, with a RMS in the solution of about 0.25 A. 
The spectra were flux calibrated with observations of the stan- 
dard stars Feige 34 and G191B2B. 

No emission features were identified (W^ < 3A) in the 
spectrum. However, an absorption doublet, identified as Mg ii 
at z = 0.351, is detected. The first line of the doublet is at 
3776.8 A, with an equivalent width (W^) of 1.3 + 0.1 A. The 
second line is at 3786.4 A, = 1.2 + 0.1 A. The doublet 
spacing is correct for the Mgii doublet at z = 0.351 (Fig. |2] 
inset). The individual line widths are not resolved however. A 
second, low significance absorption system at 35 10 A, corre- 
sponding to Fell 2600. 17 A at z = 0.351 is also marginally 
detected (Wi = 1.6 + 0.2 A). We note that the previously re- 
ported redshift (z - 0.318) was tentative and the spectrum very 
featureless, wit h lower signal-to-n oise ratio than the spectrum 
presented here jMorris et alJl99lh . 

3. Results 

3.1. X-ray lightcurves 

The EPIC lightcurve for the first observation is best fit by 
a linear increase corresponding to a rise of ~ 4% during 
the observation; the f-test yields a probability of 99.7% (3 cr) 
for the improvement in the fit compared to a constant flux 
value. The mean 0.2-lOkeV flux for this observation was 
2.80 + 0.01 X 10^'^ergcm"^ s"'. No significant variability is 



detected during the second observation, where the flux was 
3.34 + 0.02 X lO^'^ergcm^^g-i^ ig% higher than the first. 
The spectrum hardened between the observations (from F - 
2.58 + 0.015 to r = 2.41 + 0.018), consistent with previous ob- 
servations of blazars, where an increase in the synchrotron flux 
frequently coi ncides with an increase in the synch rotron peak 
frequency fe.g.lPian et al Jl99ilWehrle et al Jl998l). 

3.2. X-ray spectra 

The spectra from each epoch were found to be consistent, al- 
lowing for variability in the normalisation and slope of the 
power-laws. Both observations have therefore been fit simulta- 
neously, allowing only the power-law component of the model 
to vary independently. All X-ray spectral models included neu- 
tral absorption fixed at the Galactic level (6.28x 10^° cm"^). It is 
clear however that there is absorption above the Gala ctic level 
in the X-ray spectra of thi s source (Fig.^ see also Stocke et alJ 
ll995llWatsonetalJl999l) ; fitting a broken power-law to model 
the downturn at the soft end of the spectrum instead of absorp- 
tion results in a very flat low energy slope (F = 0.75) and a fit 
considerably worse than that obtained with neutral absorption 
atz = 0.351. 

In each observation, fitting a power-law with neutral ab- 
sorption at z = 0.351 to the combined data gives an unac- 
ceptable fit statistic ix^ = 1697 for 1476 degrees of free- 
dom, null hypothesis probability (NH) = 5 x 10"^). The data 
deviate significantly from the model at low energies, imply- 
ing that the absorption is different from that modelled or that 
the absorber is correct but that there is a difi'erent continuum 
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Table 1. Best-fit parameters and 68% confidence limits for var- 
ious fits to the combined EPIC data. The basic model used was 
a power-law with neutral absorption fixed at the Galactic level 
(6.28 X lO^'^cm"^). Added to this was either i) a neutral, solar 
abundance absorber at z = 0.351, ii) a neutral, solar abundance 
absorber at variable redshift, iii) an ionised, solar abundance 
absorber at variable redshift or iv) a neutral variable abundance 
absorber at z = 0.351. Columns 1-6 are the fitted model, ab- 
sorbing column density, redshift, abundance and ionisation pa- 
rameter, and the;^^ statistic over the number of degrees of free- 
dom for the fit. 



Model 


Nh 

(lO^o cm-2) 


z 


Z 


£ 


x' 

d.o.f. 


i) 


14.5 


0.351' 


1.0^ 


0.0' 


1696.9 




(14.0-14.9) 








1476 


ii) 


83 


2.1 


1.0^ 


0.0^ 


1587.7 




(70-97) 


(1.90-2.30) 






1475 


iii) 


89 


2.1 


1.0' 


0.003 


1586.1 




(67-74) 


(1.9-2.2) 




(0-0.011) 


1474 


iv) 


92 


0.35 1-' 


0.04 


0.0^ 


1589.4 




(81-103) 




(0.03-0.07) 




1475 



f Fixed 
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Redshift 

Fig. 3. Comparison of MOS-pn differences in the fit param- 
eters, redshift and gas absorbing column density. Confidence 
contours are 1 cr (solid line), 2 cr (dashes) and 3 cr (dots) for fits 
to the EPIC-pn (black) and combined EPIC-MOS (grey) spec- 
tra from both observations. These parameters are clearly not 
significantly affected by MOS-pn calibration differences. Inset. 
Confidence contours for the simultaneous fit of data from all 
EPIC cameras. 



below ~ 1 keV. Such a difference in continuum at soft ener- 
gies is often observed in Seyfert galaxies as a "soft excess", 
empiricall y modelled with one or more black-body compo- 
nents (e.g. 'Page et al."2002). Adding a black-body component 
to the power-law with absorption at z = 0.351 improved the 
fit significantly Or^/DoF= 1584/1474, NH = 0.02). The best- 
fit black-body temperature was 0. 1 19;^[J [Jyg keV with luminosity 
~ 8 X lO'*^ ergcm"^ s"' at z = 0.351. This model fits the data 
as well as models with variable redshift or variable abundances 
(see Table0, however the need to invoke an emission compo- 
nent in order to explain an essentially smooth absorbed contin- 
uum and the fact that a soft excess has never been reported in a 
BL Lac render it a less favoured alternative. 

In order to test the nature of the absorber it was assumed 
that the spectrum was entirely synchrotron dominated. The 
redshift, ionisation state and metallicity of the absorber were 
tested in four models: i) a neutral, solar abundance absorber 
at z = 0.351, ii) a neutral, solar abundance absorber at vari- 
able redshift, iii) an ionised, solar abundance absorber at vari- 
able redshift or iv) a neutral, variable abundance absorber at 
z - 0.351. The results of these four fits to both datasets are 
presented in Table[2 

Allowing the redshift to vary (beyond z - 0.351) signifi- 
cantly improved the fit Cr^/d.o.f. = 1587.7/1475, NH = 0.02, 
f-test probability = 4 x 10"^-^), giving a best-fit redshift of 
z - 2.1 + 0.2, though there is another comparable minimum 
in the space near z- 1.35 (Fig. |3|l. The contraint on the red- 
shift arises principally from the non-detection of a strong neu- 
tral oxygen absorption edge, placing the redshift above one. 
In the case where the absorber is ionised, there is no signif- 
icant improvement in the fit over the neutral case (x^jd.o.f - 
1586.1/1474, NH - 0.02, f-test probability compared to the 




0.15 0.2 
Relative Metal Abundance (Solar) 

Fig. 4. Confidence contours (1,2 and 3 cr) for the best-fit metal 
abundances (relative to the Solar values) and gas absorbing col- 
umn density to the EPIC spectra. EPIC-pn data are plotted in 
black, MOS in grey. The contours overlap only at greater than 
the 2 cr level. However, for both datasets the metal abundance is 
certainly less than half the solar value. Inset. Confidence con- 
tours for the simultaneous fit of data from all EPIC cameras. 



neutral absorber = 0.78), and the ionisation parameter is not 
well constrained (Table However, the lack of detection of 
strong O VII or O viii edges also places the redshift well beyond 
z - 0.351 for this model. 

Since it is the non-detection of these edges that constrains 
the redshift of the absorber in the latter cases, it is not sur- 
prising that a low abundance absorber with a redshift fixed at 
z = 0.35 1 is statistically as good a fit as the higher-redshift so- 
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lar abundance absorbers (x^/d.o.f^ 1589.4/1475, NH = 0.02). 
Allowing the redshift to vary in this case did not significantly 
improve the fit (x^/d.o.f - 1588.9/1474, f-test probability of 
0.50, compared to the model with redshift fixed at z = 0.351). 
There is some difference in the best-fit metal abundances de- 
termined from the MOS and pn datasets (see Fig. 0} due to 
calibration differences between the instruments; in particular, 
small systematic residuals may be fit as absorption edges in 
the MOS spectra. In any case, the abundance of the absorber 
must certainly be below half the solar value and is probably 
much lower (at the redshift inferred from the optical spectrum, 
z = 0.351). 

4. Discussion and conclusions 

Observations with ASCA indicated the peak of the synchrotron 
component in MS 0205.7+3509 to be betw een the UV and sof t 
X-rays at the time of those observations (Watson et al. 1999), 
consistent with results from the first XMM-Newton observation. 
The harder spectrum and higher flux during the second obser- 
vation imply a shift in the synchrotron peak frequency to higher 
energies. 

The featureless nature of the absorbed spectrum implies ei- 
ther that the absorber is metal poor or that the absorbing mate- 
rial is at a redshift much higher than the z - 0.35 1 Mg ii absorp- 
tion system detected in the optical spectrum. In either case the 
absorbing gas column density is similarly high (~ 10^^ cm"^), 
implying an optical extinction that is not observed (Ay - 5 
at z - where the gas t o dust ratio is similar to the Galactic 
value, Bohlin et al. 1978), a fact that supports the low-redshift, 
low-metallicity model for the absorber. 

It is possible that the absorber is associated with the BL Lac 
host galaxy. But the lack of gas and dust in elliptical galaxies 
jKnaPD 1999) . the typical hosts of BL Lacs, and the proxim- 
ity of the companion galaxy (2.3 ", corresponding to an appar- 
ent linear distance of ~ 11 kpc, Fal omo et alJll997 ^ strongly 
suggest a link between the compani on a nd the absorb ing gas, 
as proposed bv lFalomo et alJ ( 1 19971) and lWatson et all (|l999). 
Furthermore, this is, as far as we are aware, the highest column 
density so far observed in the spectrum of a BL Lac object (af- 
ter PKS 1413+135 which appears to have large absorption in 
its host galaxy, but the host in that case is an edge-on spiral, 

IPerlman et al 2002) . implying that the absorption seen here is 

not associated with the BL Lac host. In this case the z - 0.351 
Mg II belongs to the X-ray absorbing gas which must be poor 
in dust and metals, implying that the BL Lac is illuminating a 
fairly pristine gas cloud that is associated with the galaxy close 
to the line of sight at redshift z - 0.351. 
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